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ABSTRACT 

We have carried out 13 CO(J=2-l) observations of the active star-forming 
region N159 West in the LMC with ALMA. We have found that the CO distri¬ 
bution at a sub-pc scale is highly elongated with a small width. These elongated 
clouds called “filaments” show straight or curved distributions with a typical 
width of 0.5-1.0 pc and a length of 5-10 pc. All the known infrared YSOs are 
located toward the filaments. We have found broad CO wings of two molecu¬ 
lar outflows toward young high-mass stars in N159W-N and N159W-S, whose 
dynamical timescale is ~10 4 yrs. This is the first discovery of protostellar out¬ 
flow in external galaxies. For N159W-S which is located toward an intersection 
of two filaments we set up a hypothesis that the two filaments collided with 
each other ~10 5 yrs ago and triggered formation of the high-mass star having 
~37 Mq. The colliding clouds show significant enhancement in linewidth in the 
intersection, suggesting excitation of turbulence in the shocked interface layer 
between them as is consisten t with the magneto-hydro-dynamical numerical sim¬ 


ulations fllnoue fe Fukui 


2013! ). This turbulence increases the mass accretion rate 


to ~ 6 x 10 4 M 0 yr 4 , which is required to overcome the stellar feedback to 
form the high-mass star. 


Subject headings: ISM: clouds — ISM: molecules - ISM: kinematics and dynamics 


stars: formation 
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1. Introduction 

High-mass stars are influential in galactic evolution by dynamically affecting and 
ionizing the interstellar medium, and also by chemically enriching heavy elements via 
supernova explosions. It is of fundamental importance to understand the physical processes 
in the evolution of molecular clouds where high-mass stars are forming. There have 


e-g-, 


Zinnecker & Yorke 

— 

O 

O 

CM 

Tan et al. 

2014) 


2 0141) . In spite of these works we have not yet 


understood how high-mass star formation takes place. One of the promising candidates 
where young high-mass stars are forming is th e very dense and massive cores such as 


infrared dark clouds in the Milky Way flPeretto et al 


20131 ). Another possible candidate 


is the compressed layer formed in cloud-cloud collisions. Observations of a few super star 


clusters and smaller H II regions in the Milky Way have sh own s i gns o f tri ggere c 


formation 


of hip 

'h-mass stars in the collision-compressed layers (e.g., 

Furukawa et al. 

2009; 

Torii et al. 

2011, 

2015; 

Fukui et al. 

2014, 

2015 

). Magneto-live 

ro-dynamical numerical simulations of 


two colliding molecular clouds by 


Inoue fe Fukui 


2013)) have shown that turbulence is 


excited and the magnetic field is amplified in the collision-shocked layer between the clouds. 
These turbulence and magnetic field increase the mass accretion rate, favoring high-mass 
star formation. 

The difficulty in studying young high-mass stars lies in the considerably small number 
of young high-mass stars as compared with low-mass stars in the solar vicinity; this is in part 
due to the lower frequency of high-mass stars and the heavy sightline contamination in the 
Galactic disk. ALMA is now opening a new possibility to explore high-mass star formation 
in external galaxies by its unprecedented sensitivity and resolution, having a potential to 


revolutionize our view of high-mass star formation. The 


(LMC and SMC), at a distance of 50kpc flSchaeferl 120081 ) and 61kpc flSzewczvk et al. 


Large and Small Mag ellanic Clouc 


2009 ). 
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are actively forming high-mass stars. The LMC is an ideal l aboratory to see the evolution 


of stars and clouds thanks to the non-obscured face-on view (Subramanian fe Subramanian 


2010 

) of all the GMCs in a single galaxy (for a review 

Fukui & Kawamura 

2010) 


12 


CQ ('./=! 0 ) survey for GMCs with NANTEN 4m telescope f Fukui et al 


2001 


Yamaguchi et al 


2001 


Fukui et al 


19oJ : 


Mizuno 


20081) provided a sample of nearly 300 GMCs 


at 40 pc resolution and led to an evolutionary scheme from starless GMCs (Typ e 1) to 


active star-forming GM Cs (Type III) over a timescale of 20 Myrs flFukui et al 


Kawamura et al. 


1999 


20091 ). Aiming at revealing the hner-scale details of the molecular gas 


in the LMC, we have commenced systematic CO observations by using ALMA at sub-pc 
resolution. 

Among the nearly 300 GMCs over the LMC obtained with NANTEN, N159 is the 


brightest one with H II regions. Infrared studies h ave revea. 


active in star formation. 


Mizuno et al. 


e d nea rl y twenty yo ung high- mass 


stars in N159 with Spitzer and Hersche 

(Chen et al. 

2010, 

Wone et al. 

2011, 

Carlson et al. 

2012 

and references therein; 

Seale et al. 

2014 

), where two clumps N159 East and West are 


( 20101 ) showed that the CO J=4-3/J = l-0 ratio 


shows enhancement toward the molecular peak without a well-developed H II region in N159 
West (N159W). This high excitation condition suggests that N159W is possibly on the 
verge of high-mass star formation, and thus the initial condition of high-mass star formation 
may still hold. The preceding observations with Australia Telescope Array (ATCA), while 
low in resolution (HPBW~ 6"), presented some hint of small-scale clumps and filaments in 


N159W(|Seale et al. 


20121 ). N159W is therefore the most suitable target for the purpose of 


witnessing the oneset of high-mass star formation. 

We present the first results of the ALMA observations of N159W in this Letter mainly 
based on the 13 CO( J=2-l) data. 
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2. Observations 

We carried out ALMA Cycle 1 Band 3 (86-116 GHz) and Band 6 (211-275 GHz) 
observations toward N159W both with the main array 12m antennas and the Atacama 
Compact Array (ACA) 7m antennas. The observations centered at (aj 2000 .cn <5j2ooo.o) = 
(5 h 39 m 35®34, -69°45'33 / .'2), were carried out between October 2013 to May 2014. The target 
molecular lines were 13 CO(J=1-0), C 18 O(J=l-0), CS(J=2-1), 12 CO(J=2-l), 13 CO(J=2-l) 
and C 18 0(J=2-1) with a bandwidth of 58.6 MHz (15.3 kHz x 3840 channels). We used 
a spectral window for the observations of the continuum emission among the four with 
a bandwidth of 1875.0 MHz (488.3 kHz x 3840 channels). The radio recombination lines 
of H30a and H40a were also included in the windows. The projected baseline length of 
the 12m array ranges from 16 m to 395 m. The ACA covers 9 m to 37 m baselines. The 
calibration of the complex gains was carried out through observations of seven quasars, 
phase calibration on four quasars, and flux calibration on five solar system objects. For the 
flux calibration of the solar system objects, we used the Butler-JPL-Horizons 2012 model 
(https://science.nrao.edu/facilities/alma/aboutALMA/Technology/ALMA_Memo_Series/ 
alma594/abs594). The data were reduced using the Common Astronomy Software 
Application (CASA) package (http://casa.nrao.edu), and visibility imaged. We used the 
natural weighting for both the Band 3 and Band 6 data, providing synthesized beam sizes 
of ~2'.'5 x 1'.'8 (0.6 pc x 0.4 pc at 50kpc) and ~1"3 x O'/8 (0.3 x 0.2 pc), respectively. The 
rms noises of molecular lines of Band 3 and band 6 are ~40 mJy beam -1 and ~20 mJy 
beam -1 , respectively, in emission-free channels. The comparison of the cloud mass derived 
from the ALMA observation with that of a single dish observation described in Section 13.11 
suggests that the missing flux of the present ALMA observation is not significant. 




3. Results 


3.1. A complex filamentary structure 


Figure [T] shows the 13 CO velocity integrated intensity image of the ,7=2-1 transition. 
The distribution of the 13 CO emission is highly filamentary. The filaments, having often 
straight or curved distribution, have a typical length of 5-10 pc and a width of 0.5-1.0 pc 
defined as a full-width of the emission area at the 3<x level of the intensity integrated over 
a range of 234 to 240kms _1 , which may be anal ogous to the dominance of filaments in 


the interstellar medium of the solar vicinity (e.g., 


Andre et al. 


2013 


Molinari et al. 


2Q10J), 


possibly suggesting that filaments are ubiquitous in other galaxies as well. More details of 
the filaments will be published separately. The most active star formation is found in two 
regions as denoted by N159W-N and N159W-S in Figure [TJ both of which are associated 
with enhanced 13 CO emission. 

The cloud mass is estimated from the 12 CO(.7=2-1) intensity by assuming a conversion 
factor from th e 12 CO(J=1-0) intensity to the column density of X(CO)=7 x 10 2 °cm -2 


fjFukui et al 


20081 ) and the typical 12 CO( J =2-l)/ 12 C O (.7=1 -0) ra tio toward Hu regions 


of 0.85 (the ratio in the Orion-KL region of 


Nishimura et al. 


20151 ). We also assumed the 


absorption coefficient per unit dust mass at 1.2 mm and the dust-to-gas mass ratio to be 
0.77cm 2 g _1 and 3.5 x 10~ 3 , respectively to derive the gas mass from the dust emission 
(Herrera et al. 2013). In total, the filaments have molecular mass of 2.4 x 10 5 M 0 in 
N159 W corresponding to 35 % of the total mass w hich is estimated by the lower resolution 


study flMinamidani et al. 


2008 


Mizuno et al. 


2010). We define the N159W-N and N159W-S 


clumps at the 5 a level of Band 6 continuum (white contours in Figure 1), and the masses of 
these clumps are estimated to be 2.9 x 10 4 A7 0 and 4.1 x 1O 3 M 0 , respectively, by assuming 
the dust temperature 20 K. Their masses derived from the CO emission are 1.5 x 10 4 M 0 
and 4.2 x 1O 3 M 0 , respectively, as is consistent with the dust-emission estimate. 
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3.2. Outflows 


We have discovered two molecular outflows having velocity span of 10-20 km s _1 
in 12 CO(J=2-l). Figure |2] shows the distribution of the outflow wings. One of them 
corresponds to N159W-N and the other N159W-S. The N159W-S outflow has red-shifted 
and blue-shifted lobes which show offsets of 0.1-0.15 pc from the peak of the continuum 
emission. The outflow axis is along the east-west direction. The N159W-N outflow has only 
the blue-shifted lobe which shows an offset of 0.2 pc from the 13 CO peak. It is possible 
that the complicated gas distribution around N159W-N may mask the possible red lobe. 
The size of the red-shifted and blue-shifted lobes is less than the beam size 0.2 pc x 0.3 pc 
and the upper-limit timescale of the outflow is roughly estimated to be 10 4 yrs. This is the 
first discovery of extragalactic outflows associated with a single protostar. The positions 
of outflows in N159W-N and N159W-S coincide with YSOs identified based on the Spitzer 
data: 053937.56-6 94525.4 (hereafte r YSO -N; Chen et al. 2 010) and 053941.89-694612.0 


(hereafter YSO-S; 


Chen et al. 


2010 : P2 in 


Jones et al 


2005|), respectively. 


3.3. YSO characteristics 


Two YSOs associated with outflows, YSO-N and YSO-S, h ave been stu die d ext ensively 


at near- to far-infr ared, submillimeter, and radio wavelengt hs (ICarlson et al. 


references t 


(12006 


lerem; 


Seale et al. 


2014 : 


Indebetouw et al 


20041 ). Using the 


2012 


and 


Robitaille et al. 


20071 ) YSO model grid and spectral energy distribution (SED) fitter, we model 
all the available data including the Spitzer and Herschel fluxe s (1.2-5 00 urn), as well as 


photometry we extracted from Spitzer/IBS spectra (5-37 /mi; 


Seale et al. 


20091) . and the 


fit of the SEDs indicates that both YSO-N and YSO-S are Stage 0/1 YSOs. The mass 
and luminosity are estimated to be 31 ± 8M 0 and (1.4 ± 0.4) x 10 5 L & for YSO-N, and 
37 ± 2 M 0 and (2.0 ± 0.3) x 10 5 L & for YSO-S. These results are consistent with those from 
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Chen et al. (2010), who also used Robitaille fitter but without the Herschel constraints. 


The d ynam ical ages o: 


fitter flRobitaille et ah 


the two outflows are consistent with the ages output from the SED 
200611 . ~10 4 yrs. 


According to 3 cm radio continuum measu rements YSO-N is determined to be an 05.5V 


star, whereas YSO-S is not detected fllndebetouw et al. 


20041 . suggesting that YSO-S is in 


an earlier evolutionary state than YSO-N. T his is consistent wi th a n on-detection of the 


2.113/ini and with the weak Bry in YSO-S (ITestor et al. 


2006 b The 


Testor et al. 


(2006) 


He 


near-IR VLT data revealed that YSO-S consists of at least two sources, whereas their 
detailed physical properties and relation with the mid-/far-infrared source are yet unknown. 


3.4. Filamentary collision in N159W-S 

In Figured] N159W-N shows complicated 13 CO distribution, whereas the source 
N159W-S shows relatively simple 13 CO morphology. The 13 CO distribution in N159W-N 
consists of several filaments which are elongated generally in the direction from the 
northeast to southwest, and N159W-S is located at the tip of a V-shaped distribution of two 
filaments. We shall focus on N159W-S in the following to describe the filament distribution 
and the high-mass young star, because the simple morphology allows us to understand the 
physical process unambiguously. 

Figure [3] shows the two filaments toward N159W-S (Figure EKa) the whole velocity 
range, Figure [3](b) red-shifted filament, and Figure [3](c) blue-shifted filament, respectively). 
The two filaments overlap toward N159W-S, where the 13 CO intensity and linewidth are 
significantly enhanced. Figures |3](d-h) show position-velocity diagrams taken along the 
two filaments. We see that the filaments have small velocity span of 3kms _1 in the north 
of N159W-S which shows significantly enhanced velocity span of 8kms _1 at the 15% 
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level of the 13 C0 peak in Figure[3](g). An HST image at near infrared indicates that the 
red-shifted filament is extended toward the south beyond N159W-S (Carlson et ah 2015 
in preparation), while no CO emission is detected there in our 12 CO or 13 CO observations 
with ALMA. We also find that the blue-shifted filament has its extension beyond N159W-S 
in 13 CO. So, although the filaments are apparently terminated toward N159W-S, they are 
actually more extended, placing N159W-S in the intersection of the two filaments. 

In N159W-S the longer red-shifted filament in the east is highly elongated and mildly 
curved, having a length of 10 pc, while the other blue-shifted filament in the west is straight 
and elongated by 5 pc. N159W-S clearly demonstrates that a high-mass YSO with bipolar 
outflow is formed toward the intersection between the two thin filaments, and the velocity 
dispersion is significantly enhanced in the intersection. 


Based on these results we set up a hypothesis that formation of N159W-S was triggered 
by the collision between the two filaments. We first describe a possible scenario for 
N159W-S and then discuss the observational constraints on the collision and high-mass 
star formation. The two crossing filaments overlapping toward N159W-S give direct 
support for the present scenario. The lower limit for the relative velocity in the collision 
is given by the velocity difference of the two filaments 2-3kms _1 . The actual collision 
velocity should be higher than 2-3 km s" 1 because of the projection effect. According 
to the m agneto- hy dro-dy namical numerical simulations of two colliding molecular flows 


by Inoue & Fukuil fl2013l h the collision-shocked layer enhances isotropic turbulence, 
independent of the direction of the collision, and the velocit y span in the sh ocked layer is 


similar to the relative collision velocity. The simulations by 


Inoue fe Fukuil (120131 ) for a 


velocity difference of 20 km s _1 allow us to scale the relative velocity to ~10 km s" 1 with 
basically the same physical process. We therefore assume the velocity span in N 159W-S, 8 
km s -1 , as the actual collision velocity. This implies the relative motion of the two filaments 
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is nearly vertical, roughly 70°, to the line of sight. 


4. Discussion on the high-mass star formation processes 


Since the rest of the filaments show no sign of velocity dispersion enhancement with 
high-mass star formation, we assume that the non-interacting filaments retain the initial 
condition prior to the collision. The line-mass, mass per unit length, in the filaments 
changes from region to region by an order of magnitude. In order to estimate the typical 
mass of the filaments associated with N159W-S clump for the following discussion, we 
pick up two segments of 1.5 pc and 1.8 pc in length and 0.7pc and 0.6 pc in full width at 
a 35 % level of the 13 CO peak for the red-shifted and blue-shifted filaments, respectively, 
as indicated in Figures 3(b) and 3(c). Below the 35% level, it is hard to estimate the 
line-mass of the individual filaments separately due to overlapping. Above this level, the 
mass sampled becomes underestimated. We estimate the total mass of these two segments 
to be 2.9 x 10 3 M 0 from 12 CO(J=2-l) for a velocity range 234-240 km s” 1 . We then 
estimate the average line-mass of these two filaments to be 8.9 x 10 2 M 0 pc' 1 . The filaments 
are not detected in the Band 6 continuum at the 3 a noise level of the molecular mass 
density 1.6 x lO 3 M 0 pc' 1 , which is higher than the above CO-based line-mass density of 
the filaments. 


This suggests that the collision took place in a timescale, ~0.5pc divided by 8kms _1 , 
i.e., ~6 x 10 4 years ago. We assume that formation of the high-mass star initiated at 
the same time. By using the stellar mass 37 M 0 , the average mass accretion timescale of 
the star formation is given as 37 M 0 /6 x 10 4 yrs ~6 x 10 -4 M 0 yr' 1 . This rate is well 
in accord with the theoretical estimate around 10' 3 M 0 yr _1 and satisfies t he criterion to 


form high-mass stars by overcoming the stellar feedback (e.g., 


Wolfire fe Cassinelli 


19861 ). 


The small outflow timescale 10 4 yrs is consistent with this picture involving rapid high-mass 
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star formation. 


The present case of N159W-S has shown that the high-mass star having 37 M 0 is 
formed in a turbulent condition created by the collisional shock. The mass of the N159W-S 
clump is estimated to be 4 x 10 3 M e toward its CO peak. There is no sign of such dense 
clumps over the rest of the filament according to our ALMA data, either in CS(J=2-1) 
data, whose line-mass detection limit is about 15OM 0 pc _1 , or in dust emission data, whose 
line-mass detection limit is about 1.6 x lO 3 M 0 pc -1 by assuming a filament width of 0.6 pc. 
This offers an interesting possibility that high-mass stars do not necessarily require dense 
cloud cores as the initial condition. Instead, high velocity colliding molecul ar flow s are 
able to efficiently collect mass into a cloud core non-gravitationally. Inoue fe Fuku il (120131) 


discuss that the mass flow in the collision can be efficiently converged into a shock-induced 
core due to the oblique shock effect, and that self-gravity is not important in the beginning 


of the high-mass star formation, while soon later, in the shoe 
will play a role to form the stellar core (see also 


r. in the shoe 

s-coll 

Vaidva et al. 

2013) 


ected core, self gravity 


In the Milky Way we see increasing observational evidence for cloud-cloud collisions 
which trigger high-mass star formation. Four super star clusters, Westerlund2, NGC3603, 

isions between two clouds 


RCW38 and DS 

3(20031179, are found to 

be formed by coll 

isions 

(Furukawa et al. 

2009; 

Oliama et al. 

2010; 

Fukui et al. 

2014 

2015) 


2015). Isolated O stars with 


H II region, M20, RCW120 etc., are also suggested to be triggered by cloud-cloud collisions 


(jTorii et al 


2011 


2015)). N159W-S is in the very early stage of star formation as indicated 


by the non-detection of ionized gas, as well as by the collision scenario and SED models. 
Therefore N159W-S is an optimal source to study filamentary collision leading to star 
formation. It has been shown that the youngest O stars are formed coevally in a duration 
of < 10 5 yrs in NGC3603 and Weste rlundl by careful measurements of stellar ages with 


HST and VLT (Kudryavtseva et al. 


2 0121 ). We have here an independent estimate of the 
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stellar age by taking the advantage of the simple cloud morphology in JN 


present time-scale estimate is consistent with that of 


Kudryavtseva et al 


159W-S, and the 


( 2012 ). 


5. Conclusions 

In this Letter we presented the 13 CO (J = 2-1) observations with ALMA on the active 
star-forming region N159 West in the LMC. We have found the first two extragalactic 
protostellar molecular outflows toward young high-mass stars, whose dynamical timescale 
is ~10 4 yrs. One of the two stars N159W-S is clearly located toward the intersection of 
two filamentary clouds. We set up a hypothesis that two filaments collided with each other 
~10 5 yrs ago and triggered the formation of the high-mass star. The results demonstrate 
the unprecedented power of ALMA to resolve extragalactic star formation. 
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by ESO, AUI/NRAO and NAOJ. This work was supported by JSPS KAKENHI grant 
numbers 22244014, 22540250, 22740127, 23403001, 24224005, 26247026; by JSPS and by the 
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Fig. 1.— Velocity-integrated intensities of 13 CO( J=2-l) toward N159W with the 12-m array 
observation is presented. White and red crosses show the YSOs (Chen et al. 2010) and the 
Band 6 continuum peaks of the ALMA observation, respectively. White contours are plotted 
at the 5 o level of the Band 6 continuum emission. 
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Fig. 2.— Distributions of high-velocity wings from two outflows: (a) and (b) for YSO-N 
and (c) and (d) for YSO-S. (a, c) Image and white contours show the Band 6 continuum 
emission. Contours are plotted every 0.9 Jy beam” 1 . Red and blue contours show images of 
12 CO(J=2-l) integrated over the velocity ranges shown in red and blue dashed lines in (b) 
and (d). Contours are plotted every 1.1 Jy beam” 1 km s” 1 . White cross depicts the Band 6 
continuum peak, (b, d) The red and blue lines show averaged spectra of the outflow wing 
over the regions inside the red and blue contours, respectively, in (a) and (c). 
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Fig. 3.— Spatial and velocity distributions of two filaments in N159W-S. (a) 13 CO(J=2- 
1) intensity map integrated over a velocity range 234 -240 km s” 1 . The X offset in R.A. 
(J2000.0) and Y offset in Dec. (J2000.0) are from (4 h 39 m 41!88, -69°46T2"2). (b, c) The 
13 CO(J=2-l) map integrated over a velocity range between the dashed lines of (d) and (e), 
respectively. Magenta and cyan lines show the regions used for mass estimate in Section 
4 for the velocity range same with that in (a). The arrows depict the lengths and widths 
of the regions. The black/white contours plotted in (a)-(c) indicate the 5cr level of the 
Band 6 continuum emission, (d, e) 13 CO( J=2-l) position-velocity diagrams along the region 
between the dashed lines in (b) and (c), respectively, where the central velocity is 237 km s” 1 . 
Contours are plotted at every 6 a from 3 a (1 a = 0.013 Jybeam” 1 pc), (f-h) Position-velocity 
diagrams of 13 CO(J=2-l) along the dashed lines in (a). 

















































